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Summary 


Methods  are  given  for  obtaining  the  ambiguity  characteristics  of 
various  types  of  high-resolution  radar  waveform,  amplitude  or  phase 
modulated,  pulsed  or  CW.  The  analysis  deals  mainly  with  random 
modulations,  to  which  many  complicated  types  of  waveform  approximate. 
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1 .  Introduction 

Since  originally  proposed  as  a  useful  description  of  the  properties 
of  a  radar  modulation  (ref.  l),  the  ambiguity  function  has  been  subject 
to  conaiderable  study  and  further  development  (refs.  3-8),  but  it  is  not 
an  easy  function  to  compute.  This  paper  sets  out  the  mathematics  upon 
which  some  simple  constructions  may  be  based  and  derives  formulae  and 
sohematio  diagrams  for  certain  primitive  modulations.  Most  of  these 
are  of  randan  or  quasi -random  type,  being  the  aimplest  forma  of  high 
resolution  signal  to  treat.  They  include ,  as  a  particular  case,  the 
commonest  of  all  modulations,  the  incoherently  phased  pulse -train, 
reached  in  section  11  of  the  paper. 

As  s  summary  of  the  definitions  from  which  ws  start,  we  have  in  the 
first  place  the  generalised  auto-correlation  function  of  any  complex 
modulation  u(t).  This  is  given  by 

(la)  A  (t  ,  f  )  =  I  p--sfc)a.(p+^t)e"2'",'fpcip 

or  tha  equivalent  formula 

OP 

(lb)  A(t,f)  •  e  |  u.(p -1- t)  dp 

—  £? 

which  is  a  complex  function  of  time  t  and  frequency  f,  though  in  simple 
casus  it  often  turns  out  to  be  purely  real.  The  function  deaoribeb  the 
characteristic  response  of  an  ideal  signal  filtering  system  designed  to 
indicate  the  time-shift  t  and  frequency-shift  f  imposed  on  a  transmitted 
radar  signal  as  a  result  of  echoing  from  a  point  target.  Its  form  is 
characteristic  of  the  radar  modulation  u(t),  and  ideally  it  would  be  sero 
at  all  points  in  the  t-f  plana  except  at  the  origin,  which  corresponds  to 
the  target  position.  Thus  a  target  at  (T,  F)  would  produce  a  response 
A(t-T,  f-P),  but  in  reality  the  diffuse  fora  of  A  produces  indications 
not  only  at  (T,  F)  but  at  nany  other  points  besides.  The  resulting 
ambiguity  of  radar  measurement  is  mathematically  unavoidable,  h<  sever 
the  received  signal  may  be  processed,  because  the  ambiguity  resides 
fundamentally  in  the  transmitted  modulation. 

The  use  of  complex  functions  is  apt  to  cause  confusion  in  our  minds, 
but  is  essential  to  the  mathematical  treatment.  Some  reminders  may 
therefore  be  not  out  of  place.  The  actual  voltage  waveform  transmitted 
ia  assumed  to  have  the  form 

Real  parV  of  u  (t)  e 

where  f  is  the  radio  frequency  and  u(t)  is  a  low-frequency  modulation, 

real  for  pure  amplitude  modulation  and  complex  for  phase  modulation. 

The  voltage  response  from  the  filtering  system  is  given  by 

Real  port-  of  A(fc,f  )  e2"^ 


before  detection.  After  detection,  it  would  become 


D( 
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whore  D  represents  the  detector  characteriatic.  When  con al daring  inter¬ 
ference  between  signals  from  separate  targets,  strict  analysis  demands 
addition  of  A* sf  shifted  appropriately  in  t  and  f  and  properly  phased  in 
relation  to  each  other.  If,  however,  the  density  of  targets  is  so 
great  that  the  number  of  (effectively  randomly  phased)  superpositions 
is  large,  as  in  tha  case  of  clutter,  the  appropriate  funotion  for 
convolution  as  pointed  out  by  G.  R.  Whitfield  (ref.  9)  ie 

|Alt,f)|l. 


Clu  ^er  apart,  this  function  is  fundamental  in  the  mathematical  theory 
for  a  cingle-point  target,  and  is  known  as  the  unnormalised  ambiguity 
funotion.  Its  importance  is  due  to  its  two  properties  (ref. 


(2) 

(3) 

When  normalized. 


|A(t,f)!a  5  Aa(o,o) 
j  A(t,f  )ll  dt  af  =  aUo.o). 


the  ambiguity  function  has  the  form 

lA(M)ta 

AHo.cO 


which  can  be  regarded  as  the  actual  degree  of  ambiguity  on  a  ao&le  from 
0  to  1 ,  between  two  signals  whioh  differ  from  each  other  through  a  time- 
shift,  t,  and  a  frequency-shift,  f.  Tha  integral  of  tha  normalized 
ambiguity  function  with  respect  to  t  and  f  is  always  equal  to  unity, 
from  (3).  Thus,  in  view  of  (2),  the  ambiguity  of  a  modulation  must  be 
spread  over  unit  area  of  the  t-f  plane  at  the  very  least.  Random 
modulation  and  other  complicated  waveforms  are  mostly  aimed  at  reducing 
unwanted  ambiguity  by  spreading  it  more  thinly  over  much  larger  areas. 


2.  Siebert*  a  Theorem 


The  definition  of  A(t,  f)  is  more  symmetrical  in  t  and  f  than 
might  appear  from  (la).  An  alternative  form  (ref.  i )  is 


ACM) 


U  e2lupt  dp 

-  40 


where  U(f)  is  the  Fourier  Transform  of  u(t).  The  two  forms  (la)  and 
(4)  when  considered  together,  show  close  symmetry  in  time  and  frequenoy. 
Further,  Siebert' s  Theorem  (ref.  3)  states  that 


(5) 


T  lA(T,F)|*eawilfT'Ffc,<LTdF  =  |A(M)l’ 

-ml 


The  ambiguity  function  is  its  own  two-dimensional  Fourier  Transform 
with  time  and  frequency  interchanged. 


3.  The  ainglj  puls# 


-i+- 


Tha  whole  of  the  present  treatment  is  restricted,  for  reasons  of 
simplicity*  to  the  Gaussian  function  in  one  shape  or  fora.  There  is  an 
lose Hate  loss  of  generality  but  a  very  great  simplification,  enabling 
one  to  draw  rulc-of-thiab  conclusions  of  quite  wide  applicability  when 
refined  calculations  are  not  required.  In  keeping  with  this  approach,  we 
define  the  simple at  radar  pulse  modulation  as 


(6) 


v (fc)  »  k.  e 


-ot* 


which  is  a  single  (physically  unreali sable!)  pulse  centred  at  t  %  0.  By 
substituting  into  (l)  and  ohoosing  the  normalising  constant  as 


we  obtain 

(7) 


k 


1/4. 


A(t,f) 


Jill2 

2cu 


and  hano»  the  (normalifted)  ambiguity  function 

(8)  lACt.f)!1  -  e’akX*“ir%f  /<v 


which  is  an  elliptioal  Gaussian  hill  in  t  and  f.  Conveniently  we  may 
define  the  "width”  of  a  Gaussian  function  as  the  product  of  standard 
deviation  with  V(2*r),  because  this  makes  the  product  of  widths  of  v(t) 
and  its  spectrum  equal  to  unity  (see  Pig.  l).  Thus,  by  definition, 
the  signal  duration  and  bandwidth  be  come 


k) 


T  -  /rrA 

^  s  J  Cl/tT 


These  two  quantities  happen,  in  this  simple  instance,  to  be  the  widths 
of  the  ambiguity  function  along  its  principal  axes.  It  may  at  first 
seen  surprising  that  no  factor  of  root  two  has  crept  in,  since  ambiguity 
is  a  squared  and  therefore  a  narrowed  quantity,  but  the  narrowing  is 
exactly  balanced  by  the  broadening  due  to  auto- correlation. 


Schematically,  we  may  represent  tbs  basic  Gaussian  ambiguity  pattern 
(6)  by  drawing  a  rectangle  as  In  Pis.  2.  We  may  conveniently  think  of 
this  as  an  "equivalent"  ambiguity  surface,  distorted  gecsstric&lly  but 
Indicating  widths  and  volumes  correctly.  The  diagram  represents  a 
brick-shaped  hill  of  unit  altitude  where  shaded  black,  sere  elsewhere, 
end  of  unit  volume  as  WT  ■  1. 

The  single  pulse  is  the  most  basic  form  of  radar  with  minimal 
resolution  in  time  and  frequenoy. 


L  *  <7  vZ'rr 


Fir 


Fir  2. 


Gnuceiar  function 

1/  =  width  used  throughout  the  paper 


fve^oencij 


r 


Schematic  repre rents ti or.  of  Gaussian  Hjfbiruity 
function  G(T,  V/)  for  a  single  radar  pulse  of 
duration  T  and  bandwidth  XI. 

For  this  rirmle  pulse,  TV/  =  1. 


This  and  nil  rub  request  Figures  are  drawn-1  to  rar>e  scale 
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h .  The  "bed  of  nails" 


One  of  the  most  useful  idealized  wuveforms  in  radar  theory  is  a  train 
of  unit  delta -functions,  which  has  the  convenient  property  of  Fourier 
transforming  into  itself  when  the  pulse  repetition  period  is  taken  to  be 
unity.  From  this  waveform  any  periodic  function  can  immediately  be 
generated  by  scaling  and  convolution.  It  has  been  shown  elsewhere  that 
such  a  waveform  gives  rise  to  a  generalised  auto-correlation  function  of 
the  form 


(10)  A  Ct  ,f  )  -  H  JL  S  (t  -n)  -  m') 

n  m 

which  has  been  described  by  Robert  Price  (ref.  6)  as  a  "bed  of  nails"  - 
sharp  side  up. 

It  is  perhaps  fortunate  that  questions  of  normalization  need  not 
concern  us  unduly,  as  they  generally  become  clear  when  physical  realism 
has  been  imposed.  The  function  (10)  would  be  particularly  troublesome 
in  this  respect,  as  the  normalization  factors  associated  with  it  are 
always  zero  or  infinite.  For  example,  the  radar  modulation  leading  to 

(10)  would  consist  of  an  infinitely  long  train  of  pulses  each  of  zero 
energy,  finite  amplitude  and  zero  duration.  And  the  associated 
normalized  ambiguity  consists  of  an  infinite  number  of  packets  of  zero 
ambiguity  at  the  integer  lattice  points  in  t  and  f.  These  physical 
absurdities  are  unombarrassing,  as  the  bed  of  nails  is  merely  a  means 
to  an  end. 

Simple  changes  of  variable  show  that  the  lattice  spacingB  can  be 
generalized,  as  in  Fig.  3,  to  become  R  in  time  and  1/^i  in  frequency, 
where  R  is  any  arbitrary  repetition  period.  The  modulation  giving 
rise  to  such  a  lattice  diagram  is  of  course  the  pulse-train 

(11)  U.(t)  =  1<  Stfc-nR') 

n 

In  the  limits  R  -  0  and  R  =  oo ,  the  ambiguity  diagram  degenerates.  In 
one  direction  all  but  a  single  row  of  nails  pass  out  of  the  pioture, 
whilst  at  the  same  time  in  the  other  direction  the  nails  crowd  together 
into  what  might  be  described  as  a  razor  edge.  This  is  also  the  limit 
of  the  single  pulse  considered  in  seotion  }  when  either  W  or  T  goes  to 
zero. 


The  modulation  (li)  is  bedevilled  by  periodio  ambiguity,  but  apart 
from  this  the  resolution  in  t  and  f  is  perfeot,  except  in  the  limits 
Just  mentioned.  When  in  these  limits  the  nails  make  contact,  the 
modulation  suddenly  becomes  an  example  of  a  radar  with  minimal 
resolution. 


5.  Convolution  of  modulations 


Let  v(t)  and  w(t)  be  two  complex  modulations, 

be  their  auto-oorrelation  functions.  Let  u(t)  be 
v  and  w,  thus 

OO 

(12)  U-(t)  =  |  v  (x.)  w  (t  -  x)  dx  . 

-aji 

We  may  wish  to  know  the  auto -correlation  function  of  u.  Substituting 
in  (la),  we  obtain 


and  let  and  A^ 
the  convolution  of 
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v*(x) 


-  f 


tij)  s^*(p  -  it  -  x)  w  (p  >  it  - e  dj  dp 


Simple  manipulations  and  changes  of  variable  enable  us  to  recast  this 
expression  in  the  fora 

<K? 

03)  AaM)  =  [  Av(x,f)  Aw(t-x,f)  ax, 

~J0 

which  is  the  convolution  of  and  with  respect  to  the  tine-dimension. 

An  exactly  similar  result  can  be  derived  in  the  frequency  dimension. 

Thus  we  have  two  useful  rules: 

Rule  1  If  two  module  cions  are  convolved  in  the  time  domain  (or  multiplied 
in  the  frequency  domain),  their  generalized  auto-correlation  functions  are 
convolved  in  the  time  coordinate. 

Rule  2  If  two  modulations  are  convolved  in  the  frequency  domain  (or 
multiplied  in  the  time  domain) ,  their  generalized  auto-correlation 
functions  are  convolved  in  the  frequenoy  coordinate. 

The  operation  of  these  two  rules  is  non- commutative. 


6.  Functional  notation 


We  now  set  up  a  functional  notation  adapted  to  the  application  of 
rules  1  and  2  and  designed  to  facilitate  the  process  of  transferring 
mathematical  formulae  to  graphical  fora.  The  notation  is  not  always 
well  adapted  to  detailed  analytical  work,  such  as  the  evaluation  of 
integrals,  and  will  be  brought  into  play  only  when  appropriate. 

(a)  Gaussian  function  of  one  variable  G-(L) 

For  reasons  which  are  not  merely  obtuse,  the  variable  itself  is  not 
stated.  G-(L)  stands  for  a  gaussian  function  of  width  L,  centred  at  the 
origin,  "width"  being  defined  by 

(14)  L  =  <T  JTx 


where  o'  is  the  standard  deviation.  The  amplitude  is  taken  to  be  unity 
in  the  middle.  Thus,  for  example,  the  waveform 


u(0  = 


-TT  (t/T)1 

e 


would  be  written 

=  CtT). 


It  is  particularly  to  be  noticed  that  the  Fourier  Transform  of  this  wave¬ 
form  u,  defined  as  usual  with  conjugate  variables  t  (seconds)  and  f 
(cycles  per  second),  is  simply 

u  =  T.  qCt'1) 

(b)  Gaussian  function  of  two  variables  G(L,  L' ) 

This  denotes  the  product  of  G(L)  in  one  variable  with  G(L’)  in 
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another.  Throughout  this  paper  the  two  variables  are  understood 
to  be  t  and  f  respectively. 

(o)  Row  of  nails  in  one  variable  l(R) 

A  train  of  unit  dalU  functions  at  regular  separations  R.  For 
example,  equation  (ll)  would  be  written  as 

c.  =  k.  UR) 

(d)  Bed  of  nails  in  two  variables  l(R,  R') 

Consists  of  unit  deltA-functions  at  the  lattice  points  (nR,  mR.') 
where  n  and  a  are  integers.  For  example,  with  variables  t  and  f 
understood,  equation  (10)  would  be  written  as 

A  :  1(1.1). 

The  more  general  case  illustrated  in  Fig.  3  is 

A  =  I  (R,R“‘) 


(e)  Convolution 


A  letter  t  or  f  under  the  etar  denotea  the  variable  of 
convolution.  For  example,  equation  (13)  would  be  written 

Aa  *  Av  *•  Aw 

t 

(f)  Arbitrary  normalising  constant  k 

It  is  convenient  to  reserve  a  latter,  such  as  k  in  this  paper, 
for  a  constant  whose  value  need  not  be  consistent  in  the  analysis. 


7.  Tha  coherent  pulse-train 

The  simplest  finite  pulse-train  for  mathematical  analysis 
consists  of  the  carrier  modulated  in  amplitude  by  two  functions 
together,  a  succession  of  gaussian  shaped  pulses  and  an  overall 
gaussian  taper  long  compared  with  the  repetition  rate.  This  may 
be  written  in  the  notation  of  the  previous  section  as 


(15) 

U.(q(w-‘)*  i(r)).q(t) 

where 

R  =  repetition  period 

T  =  duration  of  train) 

W  =  bandwidth  j  WT  »  1 

Before  wo  can  apply  either  of  rules  1  and  2,  we  require  the  generalised 
auto-eoi relation  function  of  G.  In  section  3,  it  has  been  stated  that 
the  waveform  (6)  gives  the  auto-correlation  function  (7)»  Rewriting 
in  terms  of  the  widths  defined  at  (9) ,  we  have  the  easily  proved  result 
that  the  waveform 

(16) 

u  =  k  C,  (T ) 

has  the  generalized  auto-correlation  function 
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(17) 


A  -  C,  (Ta/2,  T“'J2) 


Applying  rules  1  and  2  of  section  5,  we  thus  obtain  the  auto-correlation 
function  of  (15)  in  the  form 

(18)  A  *  (q  (vi-'JZ.'Ml)  *  I.(R,R''))>  Q(tV2,T"V2) 

fe  f 

and,  after  squaring, 

(19)  Up  *  w) *  *  c 

b  r 

which  is  illustrated  in  Figure  4,  a  well-known  diagram.  It  should 
particularly  be  noted  that  the  central  sons,  which  is  a  measure  of 
the  radar  resolution,  can  be  made  as  fine  aa  desired  by  increasing 
the  duration  T  and  bandwidth  W  of  the  radar  transmission,  assuming 
that  both  are  fully  utilised. 


8.  Complex  noiaa  modulation 

The  ambiguities  shown  in  Figure  4  may  prove  embarrassing  in  a 
praotioal  situation.  They  can  be  smoothed  out,  while  preserving  the 
central  tone,  as  illustrated  in  Figure  5,  if  the  radar  transmits  noise 
modulation  of  the  same  duration  and  bandwidth  as  before.  This 
represents  a  Gaussian  hill  of  unit  amplitude  standing  on  a  lower  wider 
Causal an  pedestal;  it  is  the  aim  in  this  aeotion  of  the  paper  to 
derive  and  ooouent  upon  this  rather  well-known  result. 

First  it  is  necessary  to  take  note  of  the  statistical  nature  of 
the  ambiguity  function.  Any  particular  high-resolution  no:  as  wave¬ 
form  possesses  a  normalized  ambiguity  function  whose  value  is  unity 
at  the  origin,  and  has  a  random  skirt,  pedestal,  or  side-lobe  pattern, 
extending  wwll  beyond  the  central  tone.  A  differen.  sample  cf  noise 
would  produoe  a  different  ambiguity  function,  the  ups  and  downs  of  the 
pedestal  differing  in  their  detailed  arrangement.  The  only  simple 
answer  ia  to  work  with  the  mean  ambiguity  function  for  an  ensemble  °f 
noise  modulations,  i.e.  the  mean  squared  auto-oorrelation  function 


(20)  |  ACkjOl^ 

This  proves  to  be  unnecessarily  restrictive  in  the  analysis,  and 
instead  of  it  we  propose  to  study  the  funotion 


(21) 


A(t,,f,)  A*( tj.fi) 


which  is  the  complex  co-varianoe  of  the  generalized  auto-oorrelation 
funotion.  By  putting  the  two  t’s  and  the  two  f's  equal,  the  mean 
ambiguity  is  immediately  obtainable,  but  (2l)  can  tall  us  more.  It 
contains  information  about  tha  structure  of  the  random  variations 
of  ambiguity  in  the  pedestal.  When  (21 )  is  zero,  the  fluctuation 
of  ambiguity  at  ,  f ^ )  can  be  regards’  as  uncorrelated  with  that 

at  (tg,  fj)»  It  may  b»  important  to  know  about  this  when  a 


T 


« - — - — - > 

A 


Fig  4.  Ambiguity  of  coherent  pulse-train 

T  =  duration  of  train  R  =  repetition  period 
W  =  bandwidth 

The  above  diagram  is  drawn  with  TW  =  15 


T 


Fig  5.  Ambiguity  Z(T,  W)  of  a  burst  of  noise 
T  =  duration  of  burst  W  =  bandwidth  of  noise 


1 
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convolution  ia  performed  upon  A,  since  simple  addition  in  "power"  is 
valid  only  in  the  absenoe  of  correlation. 

The  modulation  to  be  considered  n  this  section  will  be  of  the  form 
(22)  u,  (t)  =  v(fc).n(t] 


where  v(t)  is  Gaussian  amplitude  taper  (of.  eqn  (6)) 


and  n(t)  is  a  complex  stationary  random  function  of  time.  Its  power 
spectrum  will  be  assumed  to  be  centred  upon  zero  frequency,  to  be  of 
Gaussian  shape,  and  indeed  n(t)  will  be  assumed  to  have  all  the  usual 
properties  of  Gaussian  noise.  Thus,  when  modulating  the  radio¬ 
frequency  carrier,  n(t)  would  produce  a  typical  narrow-band  noise  wave¬ 
form  displaying  random  fluctuations  both  in  amplitude  and  phase.  ( Hsd 
we  postulated  that  n(t)  be  purely  real,  this  would  not  have  been  the 
oase,  and  the  radar  signal  would  have  been  modulated  in  amplitude  only). 

To  put  a  scale  on  tho  time-structure  of  n(t),  equivalent  to  setting 
its  bandwidth,  w»  ahall  now  assume  the  following  auto-correlation 
function  (for  ooomant,  see  Appendix  2) 

(2h)  (t)  *  n(-c)  r\*(-c  + 1)  *  . 


In  a  statistical  sense,  the  waveform  u(t)  is  now  fully  determined. 
The  constants  have  been  chosen  so  as  to  make 


(25) 


lw-11  dt 


vHt)  at 


I 


The  constants  a  and  b  control  respectively  the  duration  of  the  noise 
burst  and  the  correlation  time  of  the  noise. 


Evaluation  of  (21 ),  using  ths  definition  (lb),  prooeeds  as 
follows 

(26)  A,A* 

=  A  A 

_  [[v(%)  v  C-C+  b,)  v(<r)v(<r+  t.J)  y 

JJ 

n*(x)  n(xtt,)n  (cr)  r\*(<r  +  tx)  .  e  2'"''  ^ 


It  is  shown  in  Appendix  2  that  the  system  averaged  quadruple  produot 
la  equal  to 
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(27) 


Substituting  into  (26),  we  find  that  the  first  of  these  t0 

a  separate  produo t  of  single  integrals,  which  can  be  identified  as 


(28)  A,  A* 


=  e  * 


*(t.)  Jv(T)vfc*t.)e"2irif,T4T.  vj-(ta)JvC<r)v(<r+b1) 

This  is  the  "systematic"  contribution  to  the  mean  squared  modulus  of 
A.  reducing  to  the  square  of  the  mean,  and  it  corresponds,  as  we 
might  anticipate,  to  the  oentral  sons  of  the  ambiguity 
integrals  are  readily  evaluated,  using  (6)  and  (?)  as  the  pattern, 


Zirifi<7 


Acr 


integrals 
using  (24)  also,  we  obtain 


(29) 


kX  = 


-j  («.  +  2b)(b*+t*)  - 


Hence,  identifying  suffixes  1  and  2,  we  find 

(30)  i  A  \  =  Q  (W  \T  ') 


where  T  and  W  are  given  by 


(31) 


T  »  duration  of  noise  burst 


(32) 


V  s  bandwidth  of  noise 


i/(o  1 2l>)  /Tr 


The  two  terms  contributing  to  W  arise  because  the  amplitude 
modulation  of  the  noise  broadens  the  bandwidth  slightly.  In 
an  application  the  condition 


(33) 


b»a  (i.e.  TW»0 


would  be  likely  to  apply,  and  the  spectral  broadening  due  to 
the  Gaussian  taper  would  be  very  slight. 

The  second  term  of  (27)  also  contributes  to  (26)  and  the 
resulting  double  integral  does  not  separate  so  readily.  To 
describe  the  evaluation  as  merely  tedious  would  be  an  understatement . 
The  result,  however,  is  simple  to  write  down,  and  ia 


"14“ 


w 


Hen >3 ,  by  equating  suffixes  1  and  2,  and  remembering  to  inolude  the 
term  ^29;  already  derived,  we  obtain 


(35)  \  Al1 


-tx(«-  +  2b)  -TT1f1/o  /  ~  _afc'  -Tr,f*/(o-  +  ^) 

«  e  C 


or,  in  functional  form,  using  (31 )  and  (32), 

(36)  | Aja  ®  q  (w'\ T-')  +  (TW)''q  1t,w)  *  z(t,w)  ^ 


This  is  the  result  illustrated  in  Pigure  5.  It  shows  that  ambiguity 
is  spread  with  density  l/WT  over  a  region  corresponding  to  the  product 
of  the  duration  and  bandwidth  of  the  signal,  whilst  the  oentral  sone, 
which  is  exactly  similar  in  shape,  has  a  size  governed  by  the 
reciprocals  of  signal  bandwidth  and  duration.  In  the  limit  when 
T  ■  W  *oo,  this  radar  is  perfect. 

An  interesting  question  is  whether  (36)  can  be  said  to  be  properly 
normalised.  When  no  rand oca  element  is  present,  the  normalized  value 
of  A  and  henoe  also  of  its  square,  is  unity  at  t  =  f  =  0.  This  can 
be  assured  from  equations  (l)  by  normalizing  the  radar  modulation  so 
that 

00 

(37)  f  laLt)|*at  *  \. 

When  there  is  a  random  element  in  the  waveform,  equation  (25)  merely 
ensures 


|  A  (0,0)  |  *  1  , 

but  does  not  ensure 


I A  Co,o)  I  a  =  1. 

Indeed,  from  equation  (36),  for  example,  we  have 


(38)  \  A  C0,0)  |X  -  1  -r  (TW)*' 


The  small  second  term  (cf.  (33))  ie  due  to  the  random  pedestal  upon 
whioh  the  normalized  central  sone  stands,  and  is  the  omen  squared 
deviation  of  A  from  its  squared  mean.  However,  the  right  hand  side 
of  (38)  oan  also  be  interpreted,  from  the  theorem  of  equation  (3),  as 
the  total  volume  under  the  ambiguity  surfaoe.  It  is  now  dear  that 
the  first  term  on  the  right  of  (38)  arises  from  the  integral  of  the 
second  term  of  (36),  which  is  the  mean  volume  of  the  pedestal,  while 
the  smaller  second  term  in  (38)  is  the  volume  of  the  oentral  sone. 
Thus  neither  of  the  terms  in  (38)  oan  be  associated  uniquely  with  one 
or  other  of  the  components  of  the  ambiguity  funotion  (36).  To 
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re  turn  to  the  question  of  normalization,  it  ie  dear  that  the  fully 
normalised  function  would  be 

(is)  q  (w_' T  +  tttw  q  CT-  w> 

but  there  is  no  special  requirement  for  this  fora. 

As  for  the  fluctuations  of  the  pedestal,  inspection  of  (34)  shows 
that  these  are  effectively  de-correlated  at  separations  given  by 

(W)  -  Tjrr/l 


(W)  (-F.  -UX  =  2a(a4.2t)/(irb) 

i.e.  at  separations  whioh  are  similar  in  magnitude  and  direction  to 
the  radius  vector  of  the  limit  to  the  central  zone  of  the  diagram. 
This  would  be  expected  on  purely  intuitive  grounds,  as  the  central 
zone  describes  the  typical  resolution  of  the  matched  filter  outputs. 

It  is  interesting  to  note  that  the  central  zone  is  the  two- 
dimensional  Fourier  Transform  of  the  mean  pedestal  and  vice  versa. 
This  is  evident  from  Siebert*  s  Theorem,  and  it  must  generally  ensure 
that  the  central  zone  is  similar  in  shape  to  the  outer  boundary  of 
the  ambiguity  diagram. 


9.  Bursts  of  complex  noise 

By  applying  rules  1  and  2  of  section  5,  the  basio  ambiguity 
diagram  derived  in  the  previous  section  for  a  single  noise  burst 
can  be  generalized  for  more  complicated  signal. a.  In  this  section 

we  consider  two  cases 


(i)  a  finite  train  of  identical  complex  noise  bursts 

(ii)  a  finite  train  of  independent  complex  noise  bursts 

In  oase  (i),  it  is  neoessary  to  qualify  the  word  "identioal",  as  the 
limitation  on  the  number  of  bursts  will  be  imposed  by  means  of  a  slow 
Gaussian  amplitude  taper  extending  over  an  otherwise  infinite  train. 
Let 


(42)  W 

T 
P 
4 

For  the  present, 


bandwidth  of  the  noise 
duration  of  the  train 
duration  of  each  burst 
repetition  period  between  bursts, 
shall  assume  that 


0*3) 


FW  »  1. 


Consider  case  (i).  The  passage  from  a  single  burst  to  a  train 
of  identical  bursts  is  exactly  similar  to  the  discussion  at  equations 
(l6)-(l9)  applying  to  an  ordinary  pulse. 
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Fig  7 .  Ar.biguity  of  train,  of  independent  noise  bursts 

T  =  total  duration  of  signal  W  =  bandwidth 

P  =  duration  of  each  burst  R  =  repetition  period 
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Henoe  the  ambiguity  function  is 

(44)  (Z(P,W)  «■  I(R,R-1))  *  C(T,T-) 

t»  F 


where  2  is  the  aabiguity  function  for  a  single  burst,  i.e. 


(45) 


z(p,w)  =»  c(w-,p")  +■  (pwr'c  (P,\N) 


from  equation  (36).  The  funotion  (44)  is  shown  at  Fig.  6.  At  those 
lattioe  points  which  are  not,  to  a  crude  approximation,  fully  ambiguous, 
the  density  of  aabiguity  is  found  to  be  1/WP,  as  also  in  the  wings  of 
all  the  ambiguities. 

In  oaae  (ii),  it  is  neoessary  to  give  the  single  burst  which  we 
take  as  our  starting  point  a  duration  T,  and  then  to  multiply  this 
by  a  modulation 

(46)  CCP^ltR) 

whioh  it  an  infinitely  long  train  of  amplitude  pulses  each  of  length 
P  and  separation  1.  Henoe  the  aabiguity  funotion  is 

(47)  Z(T,W)  *  (C  (  F,  P“')  *  1(P,R")) 

f  *• 

as  illustrated  in  Figure  7.  On  this  diagram,  the  bandwidth  marked  iff 
is  slightly  greater  than  the  W  appearing  in  the  expression  (47) , 
owing  to  the  effect  of  convolution  in  frequency.  In  other  words, 
the  resultant  bandwidth  has  been  slightly  broadened  by  the  amplitude 
modulation  associated  with  P,  whioh  (unlike  that  due  to  T)  is  not 
included  in  the  V  defined  at  (32)  as  used  in  the  formula  (47). 

For  completeness,  and  for  subsequent  reference,  a  further 
aabiguity  diagram  can  be  deduced  from  the  last  expression.  Let 
us  retain  the  original  definition  of  W  given  at  (32) ,  but  now  choose 
P  so  small  as  to  make 


PW  <  1 


•  i.e.  let  the  individual  pulses  be  so  short  that  there  is  no  time  for 

the  noise  to  fluctuate  within  the  pulse.  The  signal  is  then  a  pulse- 
train  whose  bandwidth  ie  entirely  controlled  by  P.  It  is  noise- 
modulated  only  irasauoh  as  each  individual  pulse  possesses  a  random 
phase  and  amplitude,  and  the  ambiguity  diagram  will  be  found  to  be 
as  shown  at  Figure  8.  A  modulation  of  this  type  ie,  of  course, 
never  used  in  a  practical  radar  and  is  of  purely  theoretical  interest. 
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10.  Noise  amplitude  modulation 

A  radar  signal  endowed  with  random  amplitude  modulation  but  no 
phase  modulation  is  a  somewhat  artificial  idea,  but  it  possesses 
mathematical  interest  and  is  inoluded  here  for  the  sake  of  generality. 
The  amplitude  modulation  is  of  the  type  which  permits  positive  and 
negative  amplitudes  (in  fact  a  Gaussian  distribution)  and  hence 
introduces  phase-reversals  of  an  artificial  character. 

The  analysis  starts  from  equation  (26)  as  before,  though  the 
conjugate  signs  are  now  redundant  both  in  that  equation  and  in  the 
definition  (24) .  As  is  clear  from  Appendix  1 ,  an  additional  term 
must  now  be  inoluded  in  (27).  This  is  a  term 

(46)  + t -«r)4/(fc,  +  o-T) 

Putting  t(  +  X  =  the  additional  contribution  to  (26)  may  be  written 
in  full  as 


V  (x t,)  V  (x')  V  (<r)  V  (or  + k*)  . 


ij/  (x'-  or)  ^  (x'-  <r  -  b,-  e  ^1<r)  Cl-c/  Act 


and,  by  inspection,  we  see  that  this  is  exactly  the  same  as  the 
contribution  to  (26)  from  the  second  term  of  (27)  with  the  sign  of  t1 
changed.  Consequently,  an  extra  term  must  be  added  to  (35)  by 
changing  the  sign  of  t^  in  (34),  i.e. 


(49) 


1 


in  (34)  must  be  replaced  by 


(50) 


€ 


Ufc.+b,)1 


Upon  equating  the  suffixes  1  and  2,  the  resulting  ambiguity  function 
is  found  to  be 


or,  in  functional  form,  using  (31)  and  (32), 


F~i  K 
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Fig  8.  Ambiguity  of  train  of  un  correlated  pulses 

=  pulse  length  T  =  duration  of  train 

=  repetition  period 

Each  pulse  has  an  independent  random  amplitude  &  phase 


2.  (Tw)"1  (tw)'1 


noise  amulitude  modulation 


T  =  duration  of  bur rt  W  =  bandwidth 


Fir  10.  Ambiguity  Y(T,  \V)  of  a  burst  of 
random  -phase  modulation 

I'otc  the  clear  central  rtrip 
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(52)  1A|X  =  C(w",  T'1)  +  (tw)"'c;(t,w)  +  {TW)_tC(W-» 


*  *  (T,W),  s«*y 


as  illustratsd  in  Figure  9.  The  analysis  is  exact. 

That  amplitude  noise  modulation  suffers  in  comparison  with  oomplex 
noise  only  in  its  ability  to  resolve  Doppler  when  there  is  no  range 
difference  must  presumably  be  due  to  the  excessive  fades  in  an 
amplitude -modulated  signal  of  the  type  considered.  When  there  is 
a  quadrature  component,  these  fades  are,  for  the  most  part,  filled  in. 

Mathematically,  the  result  is  interesting  for  two  reasons.  If 
the  parameter  b  is  set  equal  to  sero,  t*..s  noise  will  have  no  time  to 
fluctuate  before  the  burst  is  over.  Hence  the  signal  will  oonsist 
of  a  simple  Gaua sian-ehaped  pulse  of  Gaussian  randon  amplitude,  an 
extreme  form  of  non-ergodic  modulation.  From  (31 )  and  (32)  we  see 
that  WT  is  then  equal  to  unity  as  we  should  expect.  However,  we 
find  from  (52) 


(53) 

(A)2 

-  C(t,t-') 

b-  0 

(5k) 

*  3q(t,r') 

b  =  o 

By  contrast,  for  a 
(52)  wa  have 

complex  noise 

-burst  whioh  lacks  the 

third  term  in 

(55) 

uTl 

*  2  C(T,T-') 

b  *  o 

These  results  are  due  to  the  fact  that,  for  a  real  Gaussian  random 
number  n,  we  have  the  well-known  property 

(56)  Av  (n4-)  =  3[Av(n*)]a 


whilst  for  a  complex  random  Gaussian  number,  we  have 


(57) 


Av  |nvt  =  AvO*^1)1] 

=  Av  (xv)  +  2  Av(xl)Av(y1)  +  Av(^) 

-  (3  +  1  +  3)  [AvCx1)]* 

*  2  [Av  [x*+  y1)]1 

-  2  [Av  |n1)]2 


whioh  acts  as  a  crude  check  on  ths  results  so  far  obtained. 


The  second  point  of  interest  lies  in  comparing  the  ambiguity 
function  for  amplitude  modulation  with  that  for  pure  phase  modulation, 
to  which  we  now  pass. 
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11 .  Random  phase  Modulation 

Wo  now  assume  a  modulation  of  tha  fora  given  at  equation  (22),  with 
n(t)  re-defined  to  ba 

(58)  n(0  - 

whore  z(t)  ia  a  purely  real  Gaussian  random  function.  This  will  be  a 
noiee -burst,  aaooth  and  gaussian-shaped  in  amplitude  but  randomly 
varying  in  phase .  The  desired  oo-varianoe  of  auto-correlation  (26)  will 
then  entail  evaluation  of  the  system,  average  of  the  form 


(59) 


M  ■  exf  ( i-x, ix3-y  Ctcy.) 


in  place  of  the  barred  product  in  (26).  Fortunately  this  is  directly 
obtainable  from  the  aooent  generating  function  defined  at  ( V-2)  in 
Appendix  1 . 


Putting  y1  »  y^  »  1  and  y2  *  y^  ■  -1 ,  we  have 


(60) 


M. 


where 

(61)  Q 


2x,'Xx  -  2*,x,  +  2x,X|, 


+  x*  4  2  x,x, 


2xtx4 


+• 


- 

*  K 


Here  we  have  used,  the  shorthand  notation 

-  *00 

e  *(*«-*•  t.) 

=  » c<0 

-  x  (cr  +  tx) 

Now  denoting  the  auto-oorrelation  function  of  x  by  o,  we  have  from  (6l) 
and  (62) 

(63)  M«  e  ”2c^  +  c(tl)+c(tt)+c(T-<r)+c(-e-<r+t,-tl)-c(-e-<»--fc,) 

-cix-s'  +  t,) 


To  obtain  the  ambiguity  funotion,  it  is  neoessary  to  make  approximations , 
dividing  the  t-f  plane  into  two  regions  now  to  be  considered  in  turn. 
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(i)  t1  and  t2  small 

Whan  tha  t'a  are  sero,  the  exponent  in  (63)  rani shoe  entirely,  we 
have  M  ■  1 ,  and  the  whole  prohlea  becomes  identioal  to  an  unmodulated 
buret,  i.e.  to  the  aimple  pulsa  considered  in  eection  3»  This  provides 
the  due  for  approximation  whan  tha  t'e  are  small  but  not  neoessarily 
soto ,  whioh  is  to  expand  tha  functions  appearing  in  (63)  in  tha  fora  of 
power  series  and  sea  what  happens.  Auto-correlation  functions  are 
always  even  and  may  be  expanded  in  tha  form 

(61*.)  C  (t)  *  a-,  -  o.tfc  +  --- 


More  appropriately  we  may  write 


(t)  »  c(0).  (t  .  .) 


where  B,  which  will  not  really  enter  the  present  discussion,  is  inter¬ 
preted  by  the  well-known  relation 


B2  ■  varianoe  of  power  spectrum  of  x(t)  in  ops* 


Expansion  of  *n  the  terms  in  (63)  leads  to  so  much  cancellation  that 
all  terms  in  x  and  a  vanish,  and  we  are  left  with 


M  *  e 


c(t,)4-  c(Jct) -Icio) 


which  oomes  outside  the  integrals  in  (26).  To  interpret  this  result, 
we  must  relate  the  auto-correlation  function  o,  which  applies  to  the 
random  phase  angle,  to  tha  auto-oorrelation  function  for  the  actual 
noise  modulation.  At  (24)  we  assumed 


i|/(t)  *  n  (x.)  -+  fc)  «  e" 


Substituting  from  (58),  we  now  find 


^(t)  *  eyp  (-ix(-c)  -r  ix.Ct+t)) 


and  this  con  be  evaluated  in  a  similar  manner  to  (59).  Putting 
y^  =  y^  =  0  in  (A-2) ,  we  see  that 

.  ,.s  c(t)-c(o) 


and  hence  (67)  becomes 


tha  last  equation  coming  from  (68).  Identifying  suffixes  1  and  2, 
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thls  expression  can  be  introduced  aa  a  multiplying  factor  to  the 
reault  for  an  unmodulated  burst  given  at  (8).  So,  finally, 


lai*  = 


(72)  »  CtvT'.T"') 

exaotly  a a  at  (30) .  The  aignifi canoe  of  thia  result  i»  considerable, 
ainoe  the  previous  result  for  complex  noise  at  (30)  was  only  a  part  of 
that  answer.  The  ordinary  complex  noise  burst  has,  in  addition,  a 
pedestal  term  which  is,  in  the  present  ease,  entirely  absent.  There 
is,  in  faot,  a  clear  atrip  along  the  f-axis  outside  the  oentral  zone 
given  by  (72).  This  is  fully  in  accordance  with  simple  engineering 
commonsense,  sinoe  the  random  phase  modulation  does  nothing  to  impede 
a  good  Doppler  determination  provided  that  there  is  no  unknown  time- 
delay  to  produce  a  statistical  muddle  from  temporal  auto-correlation. 
If  the  range  of  the  target  is  known,  the  randomness  of  the  phase 
modulation  is  entirely  irrelevant  and  eompenaatable . 

(ii)  t.  and  large 


When  the  t's  are  large  compared  with  the  width  of  c(t),  an  entirely 
different  method  of  evaluating  (63)  is  necessary,  because  the  power 
series  expansions  of  sub-section  (i)  cease  to  be  valid.  The  trick  is 
to  regard  the  exponent  in  (63)  as  a  function  of  x  and  <r  ,  or  more 
simply  as  a  function  of  x  - <r ,  and  consider  the  effect  of  integration 
with  respect  to  -c  and  u  .  The  exponent  is  made  of  the  terms 


(a)  -2c(0) 

(b)  +  C  It,)  +  oCtJ 

(c)  -  c  (x-<r+  t,)  -  c-ix-cr -tj.) 

(d)  +  c(t-v)  t  t(T-<r  +  tl-k1) 

For  large  t,  auto-correlation  tends  to  Eero  and  the  terms  (b)  can  be 
neglected.  The  terns  (c)  produce  "dips"  in  the  integrand  at 


%-<r  =  -t,  ,  +  fc. 


but  these  occupy  little  of  the  total  area  of  integration  and  may  safely 
be  neglacted.  Finally,  then,  we  are  left  with 

(73)  ^  ,  e  c(t- -)  c(x-<r  +  -  2c(p) 


and  again  using  (70) ,  we  find 


(74)  M  *  ’Hx-cr)  -<r  +  t(  -  fct) 

This  is  exactly  the  second  of  the  terms  at  (27),  which  gave  the  result 
at  (34),  or  more  briefly  from  (36), 
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(75)  1  A)*  «  (Tw)''c(-r,W) 

This  is  the  familiar  noisy  pedestal,  but  it  is  found  only  outside  the 
central  strip  previously  discussed. 

The  complete  ambiguity  diagram  for  random  phase  modulation, 
based  on  the  results  (72)  and  (75)  is  shown  in  Figure  10  on  page  20. 

It  might  at  first  appear  that  the  above  analysis  is  independent 
of  the  standard  deviation  of  the  aotual  phase  angle  and  the  bandwidth 
of  its  temporal  fluctuation,  as  the  results  seem  only  to  involve  the 
total  bandwidth  of  the  signal  as  a  whole.  This  is,  in  reality, 
untrue.  The  methods  of  approximation  will  be  found,  upon  careful 
analysis,  to  assume  a  "reasonable"  value  for  c(0),  i.e.  an  r.m.s. 
phase  excursion  just  large  enough  to  de -correlate  the  signal  for 
large  time-delays.  This  would  imply  phase -excursions  of  the  order 
of  one  cycle,  i.e.  little  more  than  jitter.  Comparison  of  (65)  and 
(70)  and  (68)  show 

(76)  b-  2tt1B1c(.0) 

and  it  will  be  appreciated  that  b  is  proportional  to  the  square  of  the 
overall  noise  bandwidth.  The  same  value  of  b  could  be  aohieved  by 
means  of  a  large  mean  squared  phase  excursion  o(0)  and  a  correspondingly 
smaller  rate  of  fluctuation  B.  Under  these  conditions,  it  is  found 
that  the  approximation  (i)  extends  over  a  widor  central  strip  than  before, 
as  the  interface  between  the  two  regions  is  really  determined  by  B 
rather  then  «/b.  Thus  a  wider  centr*J  atrip  would  be  oloared,  but 
range  resolution  would  be  unnecessarily  worsened.  Thera  seems  little 
point  in  doing  this,  as  it  is  wasteful  of  radar  bandwidth. 

To  summarize ,  we  have  found  that,  assuming  minimum  phase  excursion 
necessary  to  de-correlate  the  radar  modulation  at  large  time-shifts, 
the  functional  formula  fcr  ambiguity  is  given  by 

(77)  CCW’.T"'),  lk!<r' 

Y(T,W)  * 

(Tw)''c  (T,w),  It|  > 


approximately.  To  a  graphical  approximation,  which  would  be  exact 
for  our  conventional  drawings,  this  could  as  well  be  written  in  the 
form 

(78)  Y(t,w)  =  q(w_l,T*')  +  (tw)'1  [c  (t,w)  -  G  (w*1,  W )] 

Comparison  with  (.36)  and  (52)  show  the  remarkable  result 

(79)  2.  Z  *  X  +  Y 

The  ambiguity  function  for  complex  noise  is  the  sum  of  the  functions 
for  amplitude  noise  and  phaae  noise  -  strange  bed-fellowa  in  most 
mathematical  problems. 
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The  result  (77),  illustrated  at  Figure  10,  lenda  itself  to 
generalisation  in  exactly  the  sane  way  as  the  "complex  noise  burst" 
given  at  (36)  and  Figure  5.  Various  types  of  pulse-train  are  obtain¬ 
able  as  described  in  section  9  merely  by  reading  Y  in  place  of  Z.  Thus 
Figure  11  shows  the  diagram  for  a  finite  succession  of  identically  phase- 
modulated  bursts.  Figure  12  shows  bursts  of  independent  phase  modulation, 
and  finally  Figure  13  is  the  diagram  for  a  conventional  incoherent  pulsed 
radar  with  no  random  phase  modulation  except  in  the  sense  that  each  pulse 
is  triggered  at  a  random  carrier  phase-angle. 


12.  Elimination  of  "range  side-lobes" 

It  goes  almost  without  saying  that  the  roles  of  time  and  frequency 
may  be  interchanged  with  impunity.  This  turns  all  diagrams  through  a 
right  angle,  and  Figure  10  thus  suggests  a  means  of  clearing  a  strip 
along  the  t-axis  free  of  random  side-lobes.  The  signal  would  consist 
of  a  spectrum  of  oscillations  phased  at  random  but  smooth  in  amplitude 
vorsus  frequency.  In  other  words,  we  might  start  with  a  conventional 
pulse  and  then  disperse  the  pulse  by  means  of  a  filter,  flat  in 
amplitude  but  with  an  elaborate  finely-structured  phase  characteristic. 
The  power  spectrum  of  the  signal  would  be  unaltered,  and  hence  the 
temporal  auto-oorrelation  function  would  remain  as  clean  as  it  was 
for  the  original  pulse. 


13.  Skew  signals 

The  ambiguity  diagrams  discussed  in  this  paper  all  possess 
rectangular  symmetry,  but  the  method  of  convolution  described  in  section 
5  can  be  powerfully  employed  to  obtain  diagrams  for  skew  signals  such 
as  linear  frequency  modulation.  One  could  take  as  a  primitive  signal, 
for  example,  the  linear  frequency  modulation 

,.\  Trirb1 


=  rt 


ft) 


(80)  ~  e 
for  which  the  instantaneous  frequency  is 

(81)  $  = 

Substitution  into  (la)  yields 

(82)  A  *  S(f- 


which  is  an  oblique  line  of  ambiguity.  Convolution  with  this 
signal  will  shear  any  ambiguity  diagram  to  any  desired  extent. 


11*.  A  Futility  Theorem 


There  is  continued  speculation  on  the  subjeot  of  ambiguity  clear¬ 
ance.  Like  slums,  ambiguity  has  a  way  of  appearing  in  one  plaoe  as 
fast  as  it  is  made  to  disappear  in  another.  That  it  must  be  conserved 
is  completely  acoepted,  but  the  thought  remains  that  ambiguity  might 
be  segregated  in  some  unwanted  part  of  the  t-f  plane  where  it  will 
ceace  to  be  a  practical  embarrassment.  We  now  put  forward  a  line 
of  reasoning,  suggested  by  a  similar  approaoh  of  Robert  Price,  which 
should  at  least  dispose  of  any  grandiose  olearanoe  schemes. 


kVl  R/(twp) 
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The  reasoning  is  based  on  Siebert's  Theorem,  which  we  have  not 
yet  applied  in  any  way.  Since  an  ambiguity  diagram  is  its  own  double 
Fourier  Transform,  we  can  apply  a  double  convoluticn  theorem,  using 
for  example  a  Gaussian  weighting  function.  Thus,  if  H  denotes  any 
ambiguity  function,  Siebert’s  Theorem  states 

(83)  H  *+  H 

where  twiddle  denotes  transposition  of  axes  and  arrow  denotes  Fourier 
Transformation  with  respect  to  t  and  f.  From  the  property  of  Gaussian 
functions,  we  have 

(810  g  (R,L)  «*-►  C(R_I,  L~') 

and  hence  by  the  convolution  theorem  applied  in  both  dimensions,  we  have 

(83)  H  *  G(R,L)  <-*  H •  C  (  L~\  R'1) 


As  an  illustration,  consider  the  ambiguity  function  shown  at  Figure  14. 
The  dotted  rectangle  represents 

(86)  C  CR  ,  R~'  ) 


The  part  of  the  diagram  inside  the  dotted  rectangle  is 

(8?)  K.GlR.R-')  -  COj",T-') 


Transposing,  transforming,  and  comparing  with  (85),  we  find 

(88)  H  *  GCR.R'1)  =  C(T,w). 


Inspection  of  the  figure  immediately  verifies  the  truth  of  this  equation, 
even  though  we  perform  the  verification  in  terms  of  rectangles  rather 
than  ellipses.  Wherever  the  dotted  rectangle  is  moved,  it  will  enclose 
the  same  amount  of  ambiguity  until  we  reach  the  edges  of  the  diagram. 

A3  a  reductio  ad  absurdum,  suppose  that  a  much  larger  area  had  been 
cleared  of  ambiguity,  (it  would  appear  at  first  that  we  already  have  a 
factor  of  almost  four  in  hand,  but  we  must  remember  that  the  diagrams 
are  not  truly  black  and  white.)  Then  it  would  be  possible  to  argue  as 
follows. 

let  G(kR,  kR  ^  be  a  rectangle,  much  larger  than  g(R,  R-1 ) ,  but 
which  still  encloses  only  the  central  pip  of  the  diagram.  Then 

(89)  H  .  G  (kR,  UR*1)  -  G  (W’  T'1) 


and  as  previously, 

(?;>)  H  *•  G  W'R-'  )  «  CL T,\M) 


i 
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Pi  r  lj.  Ai.bifn.uty  of  incoherent  -pulse-train 

T  =  duration  of  train  R  =  repetition  period 
V/  =  bnr.dv.'idth 
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From  the  result  on  the  right-hand  side,  we  ^ee  that  the  Gaussian  function 
on  the  left  can  never  (save  at  the  edges  of  the  diagram)  fit  inside  a 
clear  area  cf  H*  But  since,  in  area  covered, 

(91)  Area  cf  C  cf  q(kR,kR-'), 


this  is  contrary  to  hypothesis. 

It  follows  from  this  line  of  reasoning  that  any  ambiguity  diagram 
with  a  small  ce.  tral  pip,  i.e.  showing  high  resolution,  mu3t  extend  over 
a  reciprocal  area  in  the  t-f  plane,  and  the  clearer  the  surround  to  the 
central  zone,  the  more  uniformly  must  be  ambiguity  be  spread  over  the 
entire  plane.  No  magic  answers  are  to  be  found. 
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Appendix  1  -  Moment  generation  fo r  correlated  normal  variates 

Let  x,  ...  x  be  a  Gaussian  random  numbers  with  zero  means,  and 

12  s 

not  necessarily  independent.  Further  assume  that  all  co-variances 

(A_1)  Ki  ~ 

are  given.  Then  the  joint  distribution  of  the  s  variables  is 
completely  determined.  However,  it  is  often  unnecessary  to  know  the 
distribution  explicitly.  The  moment-generating  function  (see  ref.  2)  is 


*3  = 


:,...xs)  e^XklJk  dxs 


where  p(x.  ...  x  )  is  the  Gaussian  probability  density  function.  To 

I  s 

obtain  any  desired  moments  of  the  x's  in  p,  we  merely  differentiate  g 
with  respect  to  the  corresponding  y* s  and  then  put  all  the  y's  equal  to 
zero.  Hence,  for  example,  with  s  =  4  we  find 


t  X.XjXjX^  = 


*-  -J  u  -r  Li 


’  z  J  -Td-Jk-Jej 

and  after  performing  the  necessary  manipulation, 

(A-3)  =  *7*4.  +  +  *^.-*1*3 

The  application  required  in  section  10  of  tho  text  is  for 


*1  - 

n  (t) 

n  (.x  +  t,) 

* 

n  i<r) 

X-W  = 

n  (<r  +  fct) 

where  n  is  a  stationary  Gaussian  random  function.  Hence 


(A-^)  y-^xixixlf  =  ^(fc|)4^(fci)  +■  +  (t-<t)iHt-0>  t-tz) 

+  <P  (t'  +  z-<r)^(b x  +  <r-x) 


as  shovm  at  (27)  and  (48). 


Appendix  2  -  Complex  nolag 

Little  aeeaa  to  hare  been  published  on  the  subject  of  oomplex  random 
funotiona;  moat  treatiaea  on  noiae  strangely  fail  to  take  full  advantage 
of  the  complex  notation  which  ia  commonplace  In  aignal  analyaia.  Vo  must 
therefore  define  with  some  care  what  ia  meant  by  a  complex  normal  variate 
and  its  generalisation,  the  complex  Gaussian  random  function. 

Let  x  and  y  be  normal  variates  of  sero  mean.  Then 

\ 

X  =  pc  -h  ^ 

will  be  called  a  oomplex  normal  variate  if  and  only  if 
(A-5)  P  -  0. 


This  implies  that  x  and  y  are  uncorrelated  but  have  the  same  standard 
deviation.  The  variance  of  s  ia  defined  to  be 


(A-6)  ZZ*  =  xx  .j1  e  2  oc1  «  2  >jX 

If  s^  and  are  two  independent  oomplex  variates,  we  have 
(A-7)  VZX  *  O 

because  all  four  terms  in  the  expansion  (x^g,  y^g,  and  x2y-()hav® 

sero  means.  However,  we  shall  show  below  that  the  above  condition  ia 
true,  not  merely  for  independent  complex  variates,  but  for  any  pair  of 
variates,  whether  independent  or  not. 

The  general  multi-variate  situation  is  set  up  by  making  a  complex 
linear  transformation  of  a  set  of  independent  oomplex  numbers  s  . . .  s  . 
Let  this  be  1  8 

(A-8)  nk  « 

where  the  a* s  are  oomplex  constants.  Then  the  n's  are  correlated 
complex  numbers,  and 

nkn»  =  X  a.  x.T.  a,.t. 

*  i  i  K'  «■  i.  It  t 

Upon  averaging,  each  term  has  a  factor  like  (A-5)  or  (A-7)  whether  k  =  l 
or  not.  Henoe 


(A-9) 


nknt  *  0 


invariably,  and  we  may  take  this  aa  characteristic  of  the  oomplex  multi¬ 
variate  distribution.  If  it  holds  good  for  any  initial  set  of 
variables,  it  must  continue  to  hold  good  after  these  have  undergone 
any  number  of  linear  transformations. 
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Correlation  between  a  pair  of  complex  numbers  must  be  expressed 
by  the  complex  co-varianoe 


(aho)  n,  n  * 

and  in  practical  applications  it  is  always  the  real  part  which  is  of 
prime  importance.  This  is  due  to  the  fact  that 

(A-l i )  In,-"*!1  =  (ni  -  ) 

*  Kl*  +  |n.f-  iQ.  (o,n*). 


Thus,  for  instance,  when  two  waveforms  are  cross-correlated  to  obtain 
the  best  position  of  fit  or  minimum  squared  modulus  of  difference,  it 
is 

^  |  u.,  ux  dfc 

which  must  be  maximised,  and  the  imaginary  part  is  not  relevant.  In 
many  practical  oases,  the  complex  co-varianoe  proves  to  be  purely  real. 
For  example,  it  is  readily  shown  that  independent  complex  normal 
variates  subjected  to  a  purely  real  linear  transformation  yield  numbers 
which  have  real  co-variances.  For  let  the  a's  in  (A-8)  be  real.  Then 
we  have 


rvn 


k"i 


*  X. 


tki'z,L 


:  X-  a  z  * 

<U  t 


and  each  term  in  the  averaged  produot  is  either  real  or  sero.  It  is 
interesting  to  note  that,  when 


nin* 

is  real,  we  have  (denoting  the  real  and  imaginary  parts  by  u  and  v)  by 
hypothesis 


(AH  2) 


u2v 


but  we  also  have 


(AH  3) 
from  (A-9)« 
(ah  k) 


U.Vj. 


Thus,  when  the 


u.^ 


+  *  ° 

co-varianoe  is  real,  we  find  that 

*  uX7.  *  0 


It  follows  from  the  above  discussion  that  stationary  white  complex 
noise,  passed  through  a  real  linear  oirouit  (stereo?),  emerges  with  a 
purely  real  auto- correlation  function.  It  is  thus  sufficient  that  the 
power  spectrum  of  the  noise  should  be  synmetric&l  about  zero  frequency, 


-36- 


and  we  may  then  safely  assume,  from  (A-14)  or  indeed  from  the  physical 
insight  provided  by  the  thought  of  a  stereo  channel  for  quadrature 
components,  that  the  real  and  imaginary  parts  of  the  noise  are 
uncorrelated  from  each  other  at  whatever  instants  of  time  we  core  to 
choose. 

At  (26)  in  the  text,  it  _s  required  to  evaluate  a  fourth  order 
complex  moment  of  the  form 

3y  expanding  into  real  and  imaginary  parts,  assuming  n  to  have  a  real 
auto-correlation  function,  and  by  using  the  results  (A-3)  and  (A-14), 
it  can  be  shown  that 


(A  - 15) 

n*ri!n3n* 

*  4',24V  + 

where 

(A-16) 

4^1 

~  nwnt* 

Vnis  is  how  (27) 
extends  to  noise 
investigated. 

was  obtained.  Whether  or  not  equation  (AH5) 
whose  co-variance  is  not  real  has  not  been 

It  was  not  required  in  the  paper. 
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